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An analytical application of Niedermayer's algorithm to the Edwards-Anderson 
model: analytical results for the multicritical point on the Nishimori line 

Chiaki Yamaguchi 

Kosugichou 1-359, Nakahara-ku, Kawasaki-si, Kanagawa 211-0063, Japan 

We apply analytically Niedermayer's algorithm to the Edwards- Anderson model on random 
graphs with arbitary degree distributions. The results for the multicritical point on the Nishi- 
mori line are shown. The results are shown by applying a criterion for spin models on the 
random graphs with arbitary degree distributions. The application of Niedermayer's algorithm 
makes the size of the Fortuin-Kasteleyn cluster small and shifts the percolation threshold. The 
results for the ±J model and the Gaussian model are respectively shown. In the present ar- 
ticle, it is respectively shown for the ±J model and the Gaussian model that, by adjusting 
an introduced parameter for Niedermayer's algorithm, the percolation threshold obtained in 
the present article agrees with the location of the multicritical point. We naively estimate the 
locations of the multicritical points for the ±J model and the Gaussian model on the randam 
graphs with arbitary degree distributions. 
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1. Introduction 

Study of spin models on complex network has been 
important 1 with the development of study of complex 
network. In the present article, random graphs with ar- 
bitary degree distributions are adopted as examples of 
the complex network. In the present article, behavior of 
spins on no growing network is investigated. 

We investigate the Edwards-Anderson model 2 as an 
Ising spin-glass model. The understandings for the 
Edwards- Anderson model on random graphs and on the 
Bethe lattice are still uncompleted. l ' 3,4 In the present 
article, the ±J model and the Gaussian model in the 
Edwards-Anderson model are respectively investigated. 
The Nishimori line in the Edwards- Anderson model is a 
line in the phase diagram for the exchange interactions 
and the temperature. The internal energy, the upper 
bound of the specific heat, and so forth are exactly cal- 
culated on the Nishimori line in the Edwards- Anderson 
model. 5-9 The location of the multicritical point in the 
Edwards- Anderson model on the square lattice is con- 
jectured, and it is shown that the conjectured value is in 
good agreement with the other numerical estimates. 10 In 
the present article, the results on the Nishimori line are 
shown. 

There is a case for occuring a percolation transition 
of network that the network is devided into a lot of net- 
works by deleting nodes and/or links on a network. We 
define the percolation transition as the percolation tran- 
sition of network in the present article. There is a case for 
occuring a percolation transition of cluster that the clus- 
ter consists of bonds put between spins becomes a giant 
component. The bond is fictitious. We define the percola- 
tion transition as the percolation transition of cluster in 
the present article. In the present article, the percolation 
transition of cluster on a complex network is mentioned. 

Powerful Monte Carlo methods, using the Fortuin- 
Kasteleyn (FK) cluster in the FK random cluster 
model, 11 that are called the cluster algorithms have been 



proposed. 12-16 The Swendsen-Wang (SW) algorithm, 12 
using the FK cluster, may be one of the most popu- 
lar cluster algorithms for classical spin models. In the 
Edwards- Anderson model that has a conflict in the inter- 
actions, the percolation transition point of the FK cluster 
by the SW algorithm disagrees with the phase transi- 
tion point. 17-19 The disagreement makes the cluster al- 
gorithm inefficient. The application of Niedermayer's al- 
gorithm 20 makes the cluster small and shifts the perco- 
lation threshold. Originally, Niedermayer's algorithm is 
proposed as a Monte Carlo method for exact numerical 
estimates. The efficiency of Niedermayer's algorithm in 
the Edwards-Anderson model is shown by Liang. 21 In 
the present article, we apply Niedermayer's algorithm as 
an exact analytical method. 

In addition to the application of a gauge transforma- 
tion, 5 ' 23 our results are shown by applying a criterion 19 
for spin models on the random graphs with arbitary de- 
gree distributions. In Ref., 19 by applying the criterion 
with the gauge transformation, the percolation thresh- 
olds of the FK cluster for the Edwards- Anderson model 
on the random graphs with arbitary degree distributions 
are analytically calculated on the Nishimori line. 

In the present article, we naively estimate the loca- 
tions of the multicritical points for the ±J model and 
the Gaussian model on the random graphs with arbitary 
degree distributions. 

The present article is organised as follows. First, 
a complex network model and the Edwards- Anderson 
model are described. Niedermayer's algorithm is de- 
scribed. A criterion for percolation of cluster is described. 
The percolation threshold is shown for the ±J model. 
The percolation threshold is shown for the Gaussian 
model. The present article is summarized. 
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Fig. 1. A relation between the aspect a of the network and the 
model on the corresponding network of the aspect a. 



2. A complex network model and the Edwards- 
Anderson model 

Network consists of nodes and links connected between 
nodes. In the present article, the network models are ran- 
dom graphs with arbitary degree distributions. The net- 
work has no correlation between nodes. The node degree, 
k, is generated according to a distribution of the node 
degree, p(k). The links are randomly connected between 
nodes. 

We define a variable b(i,j) where b(i, j) gives one when 
node i and node j is connected by a link and gives zero 
when node i and node j is not connected by the link. 
The degree k(i) of node i is given by 



fc(i) = £&(*,./)• 



(1) 



The average of the node degree for links, < k >jy, is 
given by 



1 N 

< k >n= jj^Hi) 



(2) 



where < > at denotes the average over the whole network. 
N is the number of nodes. The average of the square of 
the node degree for links, < k 2 >n, is given by 



< k 2 > 



We define 



2 < k > 



N 



< k 2 > 



N 



(3) 



(4) 



where a represents an aspect of the network. 

Fig. 1 shows a relation between the aspect a of the net- 
work and the model on the corresponding network of the 
aspect a. A network is almost the complete graph when 
a is close to zero, and the model on the network is almost 
an infinite-range model. We define a coordination num- 
ber as z. p{k) = 5k, z where S is the Dirac delta function. 
The coordination number z is two when the aspect a of 
the network is one. The network consists of many cycle 
graphs when z is two. The model on the network consists 
of many chain models. In the Erdos-Rcnyi (ER) random 
graph model and in the Gilbert model, the distribution 
of node degree is the Poisson distribution. 1 Therefore, in 
the ER random graph model and in the Gilbert model, 
< k 2 > n=< k >n (< k > N +1). The ER random graph 
model is a network model that the network consists of 
the fixed number of nodes and the fixed number of links, 
and the links are randomly connected between the nodes. 
The Gilber model is a network model that link between 



nodes is connected with a given probability. 

The Hamiltonian for the Edwards- Anderson model, W, 
is given by 



N 



u 



\J2 E 

1 {j\b(i,j)= 



J ij S%Sj 



(5) 



i} 



where Si = ±1. is the strength of the exchange inter- 
action between the spin on node i and the spin on node 
j. The value of Jij is given with the distribution P(Jij). 
The ± J model and the Gaussian model are respectively 
given with the difference of P(Jij). For the ±J model, 
the distribution P i - ±J \j lj ) of Jij is given by 



P {±J \j lJ )= P 5. h] .j + (l-v)5, hj . r 



(6) 



where J > 0. p is the probability that the interaction 
is the ferromagnetic interaction (Jy = J). 1—p is the 
probability that the interaction is the antiferromagnetic 
interaction (Jy = — J). For the Gaussian model, the dis- 
tribution p( Gaussian )(J i:( ) of J^ is given by 



p(Gaussian) 



(Jij -Jo) 



V^l 2 



(7) 



The average of Jij is [Jij]n = Jo where []r denotes 
the random configuration average. The variance of 

^[Jfj]H-[J ij ? R = J 2 . 

In the present article, for calculating thermodynamic 
quantities, a gauge transformation: 

S t -> S t a t (8) 



Jij ^ J ij @i &j 

is used where <Tj = ±1. Using the gauge transformation, 
H — > H and P{Jij) — > P{JijOiOj). This transformation 
has no effect on thermodynamic quantities. 23 

3. Niedermayer's algorithm 

We define a parameter £ where the range of £ is 



< t < 1. 



(9) 



£ is adjustable. The bond is put between spins with the 
corresponding probability Pbond(Si, Sj, Jij) of the spin 
states and the strength of the exchange interaction. The 
probability Pbond(Si, Sj, J^) is 



-f'bond i Sj Jij) — 1 6 J " ^ 



(10) 



where (3 is the inverse temperature and j3 = 1/ksT. ks 
is the Boltzmann constant and T is the temperature. By 
connecting the bonds, the clusters by Niedermayer's al- 
gorithm are generated. We define an index of the cluster 
on node i as c(i). We define the probability for flipping 
the cluster as Ptu„. Pf 



flip 



IS 



Pflip 



-2/3(1-4) E 



<ij|c(;)#c(j)> 



3 % ■> S i S 7 



(ii) 



where < xy > denotes the nearest neighbor pairs con- 
nected by links, and < xy \ c(x) ^ c(y) > denotes 
the nearest neighbor pairs connected by links on con- 
dition that c(x) ^ c(y). The algorithm is the Metropolis 
method 24 when £ = 0. The algorithm is the SW algo- 
rithm 12 when £ = 1. 

The thermodynamic quantity of the bond put be- 
tween the spin on node i and the spin on node j, 



J. Phys. Soc. Jpn. Full Paper 

[< b hond (i,j) >t]r, is 

[< b hond (i,j) > T ]r = [< PbondiSuS^Jij) > t ]r (12) 

where <>t denotes the thermal average. The thermo- 
dynamic quantity of the node degree for bonds at node 

i, [< fcbond(i) >t]r, is 

[< fcbond («) >t]r = 

bond (SuS^Jij) >t]r- (13) 

0|f>W)=l} 

The thermodynamic quantity of the square of the node 
degree for bonds at node i, [< fcb on d(*) >t]r, is 

[< *bond(*) >T}R = 

[< y ' ^ ' -Pbond(<Si7 Sj, Jjj) X 

0'|6(i,j')=l} {i|6(*,/)=l} 



bond 



(Sj, S;, J,/)(l - 5j t i) + 



(Si, Sj,Jij) >t]r ■ 



(14) 



The thermodynamic quantity of the node degree for 
bonds, [< fcbond >t]r, is 



1 



[< fcbond >t]r = 2^[< fcbond(«) >T\R ■ (15) 

i 

The thermodynamic quantity of the square of the node 
degree for bonds, [< fcg ond >t]r, is 



1 N 

[< fcbond >T}R = Y^l< fc bond« >t] R ■ 



(16) 



4. A criterion for percolation of cluster 

The criterion of the percolation of cluster for spin mod- 
els on the random graphs with arbitary degree distribu- 
tions is given by 

[< fcbond >T}R > 2[< fcbond >t] R (17) 

. 19 Ineq. (17) is given by the inequality when the cluster 
is percolated. Ineq. (17) is given by the equality when the 
cluster is at the percolation transition point. It is antici- 
pated that ineq. (17) is true for sufficiently large number 
of nodes in the case that the bond does not depend on 
the size of the degree k(i). 

We define a variable for the inverse temperature (3 as 
p{/3). We set 



< P (P) < 1 . 



(18) 



We consider a case that [< 6bond(«,j) >t]r, [< 
fcbond(i) >t]r, and [< kl ond (i) > T ] R are respectively 
written in 

[< b hond (l,j) >t]r =P(P), (19) 

[< fcbond(t) >t]r = p{P)k(i), (20) 

[< fcbondW >T}R = P 2 (0)k(i)[k(i)-1] + 

p{P)k(i). (21) 
In the case, it is implied that the bias for the size of k(i) 
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does not appear in the statistical results of bonds put be- 
tween spins. Therefore, in the case that [< febond >t 
}r, [< fcbond(i) >t]r, and [< fcb ond (i) >t]r are respec- 
tively written in eqs. (19, 20, 21), the bond does not 
depend on the size of k(i). 

5. The ±J model 

The distribution P( ±J \Jij) of for the ±J model 
is, using eq. (6), given by 



P {±J \j ij ) = 



2cosh(/3 P J) 



, J tj = ±J (22) 



where (3p is given by 



a 1 i p 
Pp = — In- 



2J 1 -p 



(23) 



The thermodynamic quantity of the bond put be- 
tween the spin on node i and the spin on node j, 
[< b hond (i, j) >t] r ±J \ is, using eqs. (8, 10, 12, 22) when 

P = Pp, 

[< b hond (i,j) >t} { r J) 

{Jim} <lrn> 

E{s i} Pbond(S t , Sj, Jjj) J '™ s ' s ™ 

^{Si} 



2 N [2cosh((3 P J)} N B 

E E P^iSuSj, J tj ) e*» £< im > J lm s t s m 

{J|m}{S«} 

l _ e -2/3 P €J 

= 1 + e -2p P J ( 24 ) 

where Np is the number of all links, and Np = N < 
k > N /2. We set 

1 _ e -i&pU 



(25) 



I _|_ e -2f) P J ■ 

The thermodynamic quantity of the node degree for 
bonds at node i, [< fcbond (i) >t] r ±J \ is, using eqs. (8, 
10, 13, 22, 25) when /3 = /3 P , 

[< fcbond (i) >t] r ±J) = P i±J) Wp) k{i) . (26) 
The thermodynamic quantity of the square of the node 
degree for bonds at node i, [< fc^ondW >t ^r J ^ using 
eqs. (8, 10, 14, 22, 25) when j3 = /3 P , 

[< fcbondW >t] ( r ±J) = [ P ( - ±J \p P )] 2 k(i)[k(i)-l] + 

p {±J \l3 P )k{€). (27) 

Eqs. (24, 25, 26, 27) arc formulated as eqs. (19, 20, 21), 
respectively. Therefore, the bond does not depend on the 
size of k(i). Using eqs. (15, 16, 17, 25, 26, 27), we obtain 

(2 - o)(l - e- 2 ^« J ) > o(l + e- 2f3pJ ) . (28) 

Ineq. (28) is given by the inequality when the cluster 
is percolated. Ineq. (28) is given by the equality when 
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the cluster is at the percolation transition point. From 
eq. (28), the percolation threshold is shifted by adjusting 
£ while eq. (23) holds. Therefore, it is realized that, by 
adjusting £, the percolation threshold agrees with the 
location of the multicritical point since, by adjusting £, 
the percolation threshold is shifted along the Nishimori 
line. 

We naively estimate the location of the multicritical 
point for the ± J model on the random graphs with ar- 
bitary degree distributions. We assume that £ is the func- 
tion of a. The network is at the percolation point when a 
is one. 19 The exchange interaction is only the ferromag- 
netic interaction when a is one. Therefore, £ is one when 
a is one since the algorithm is the SW algorithm 12 when £ 
is one. The percolation transition point of the cluster by 
the SW algorithm agrees with the phase transition point 
in ferromagnetic spin models. 25 We set < k >n= N — I, 



<k 2 > 



N- 



(N-iy 



2/(N - 1), and J ->• J/VlV. 



From the settings, the model on the network becomes 
the infinit-range ±J model. Using eq. (28), we obtain 



2/3p{,7 2gp.J 

(N-2)(l-e v^)>l + e vTv 



(29) 



From numerical estimation of eq. (29), £ = O(-^p) since 
(3p has a finite value only when £ = (3(-i=). Therefore, 
we naively estimate 



(30) 



Using eqs. (28, 30), we obtain 

(2 - a)(l - e - 2 v^3p./) > a (i + e -2/3pJ) _ ( 31 ) 

Ineq. (31) gives the location of the multicritical point. Us- 
ing eqs. (23, 29, 30), the location of the multicritical point 
for the infinit-range ± J model is numerically estimated 
asp* ~ 0.804- l/(l+exp(-\/2)) and Tp ~ 1.414 ~ V2 
when J/ks = 1 is set. Ineq. (31) is obtained based on 
some assumptions. We have assumed that £ is written as 
a simple form of a, and have assumed that the location 
of the multicritical point for the infinit-range ± J model 
has a finite value. 

6. The Gaussian model 

The distribution p( Gaussian )(J ij ) of J i} for the Gaus- 
sian model is given in eq. (7). We set 

Jo 



fip = 



J 2 



(32) 



The thermodynamic quantity of the bond put be- 
tween the spin on node i and the spin on node j, 
[< &bond(*,j) > T ](^ aussian ) j i Sj us i n g eqs. (7, 8, 10, 12, 
32) when /3 = /3 P , 

r , /. ... i (Gaussian) 

[< b bond {t,j) > T \ R 

fOO fOO 



/CO rOO 
••• / I] dJ i™ n ^ (GaUSSian) (^m) x 

<lm> <lm> 

E {M iWft, Sj, Jjj)e^ ^<^> J "" s ' s " 
E{Si} e^E< Ira >^ |Sm 



2 n (2ttJ 2 ) n ^ 2 ' 

/oo />oo _ 
■ ■ ■ / ] { dJlm ^ PbondjSj, Sj, Jjj) X 
-OO J —OO w r r. l 



<lm> 



{Si} 



g- E< im> ^JT+'SP E< im > JlmSlSm 

i[ 1 + e rf (^)]- 

i{l+erf[^(l-2{)](e- J «''-'> I «'-« (33) 



where erf (x) is the error function of x. We set 



^ (Gaussian) (/3p) = ^[l + crf(^|)] 



i{l + erf[^(l-20]} e - 2 ^ 2 ^.(34) 

The thermodynamic quantity of the node degree for 
bonds at node i, [< fch nd(*) > T ]( Gaussian ) ) is, using 
eqs. (7, 8, 10, 13, 32, 34) when /3 = fi P , 

[< fcbondW > T ]< ? Ga ~) = p (G—n) {M k{{) (35) 

The thermodynamic quantity of the square of the node 
degree for bonds at node i, [< k 2 ond (i) >t]^ 
using eqs. (7, 8, 10, 14, 32, 34) when /3 = /3 P , 

7 2 /-\ - 1 (Gaussian) 
[< fc bond(«) >t\r 

[j0 (Gaussian) (/ o p)] 2 fc(z)[fc(z) _ 1] + 



(Gaussian ) 



, IS, 



(Gaussian) 



(pp)k(i) 



(36) 



Eqs. (33, 34, 35, 36) are formulated as eqs. (19, 20, 21), 
respectively. Therefore, the bond does not depend on the 
size of k(i). Using eqs. (15, 16, 17, 34, 35, 36), we obtain 

(2 - a){l + erf[^(l - 2£)]} e -2(0pJ) a «i-O < 
v2 



(2-a)[l + erf(^)]-2a. 



(37) 



Ineq. (37) is given by the inequality when the cluster 
is percolated. Ineq. (37) is given by the equality when 
the cluster is at the percolation transition point. From 
eq. (37), the percolation threshold is shifted by adjusting 
£ while eq. (32) holds. Therefore, it is realized that, by 
adjusting £, the percolation threshold agrees with the 
location of the multicritical point since, by adjusting £, 
the percolation threshold is shifted along the Nishimori 
line. 

We naively estimate the location of the multicritical 
point for the Gaussian model on the random graphs 
with arbitary degree distributions. We assume that £ 
is the function of a. The network is at the percolation 
point when a is one. 19 The exchange interaction is only 
the ferromagnetic interaction when a is one. Therefore, 
£ is one when a is one since the algorithm is the SW 
algorithm 12 when £ is one. The percolation transition 
point of the cluster by the SW algorithm agrees with 
the phase transition point in ferromagnetic spin mod- 
els. 25 We set < k > N = N - 1, < k 2 > N = (N - l) 2 , 
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a = 2/(N - 1), J -5" Jo/N, and J -> J/y/N. From 
the settings, the model on the network becomes the 
Sherrington-Kirkpatrick (SK) model. 22 Using eq. (37), 
we obtain 

(N - 2){1 + erf [^=(1 - 2£)]} e~ " < 



(AT- 2 )[l + erf(^£)]-2. 



(38) 



From numerical estimation of cq. (38), £ = 0(^=) since 
/3p has a finite value only when £ = 0(-^=). Therefore, 
we naively estimate 

? = vz. (39) 

Using eqs. (37, 39), we obtain 

(2 - o){l + erf [%^(1 - 2>/5)]} x 



V2 

e -2( / 3 P J) 2 v^(l-v^) < 

(2-a)[l+erf(^)]-2a. 



(40) 



Incq. (40) gives the location of the multicritical point. 
Using eqs. (32, 38, 39), the location of the multicritical 
point for the SK model is numerically estimated as Jq ~ 
0.886 - y^T/Z and T* P ~ 1.128 - v / 4/^ when J/k B = 1 
is set. Ineq. (40) is obtained based on some assumptions. 
We have assumed that £ is written as a simple form of a, 
and have assumed that the location of the multicritical 
point for the SK model has a finite value. 

7. Summary 

In the present article, the analitycal result of the ap- 
plication of Niedermayer's algorithm is obtained on the 
Nishimori line. It is respectively shown for the ± J model 
and the Gaussian model that, by adjusting the intro- 
duced parameter £ for Niedermayer's algorithm, the per- 
colation threshold obtained in the present article agrees 
with the location of the multicritical point. In the present 
article, we have naively estimated the locations of the 
multicritical points for the ± J model and the Gaussian 
model on the random graphs with arbitary degree dis- 



tributions. The locations of the multicritical points are 
obtained based on some assumptions. We have assumed 
that £ is written as a simple form of the aspect a of the 
network, and the locations of the multicritical points for 
the infinit models have finite values. 
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